Bifurcation, mode coupling and noise in a nonlinear multimode superconducting RF 

resonator 
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The addition of a 7-SQUID array to a thin film superconducting coplanar waveguide resonator 
introduces a controlled, tunable, nonlinear inductance that strengthens effects analogous to the Kerr 
effects in quantum optics. Exploiting the self-Kerr effect in one resonant mode, we create a sensitive 
bifurcation amplifier and then demonstrate spectroscopy of other resonant modes by simultaneously 
employing the cross-Kerr effect. We find that the sensitivity of the device approaches quantitative 
agreement with a quantum model by Dykman, but is limited by a noise that has its source(s) 
on-chip. 

PACS numbers: 85.25.-j, 74.40.-n, 42.65.Hw, 78.70.Gq 



Superconducting circuits that involve the quantisation 
of charge or flux have been the subject of great attention 
in recent years as potential building blocks for quantum 
computers. With the introduction of circuit QED [TH3], 
the controlled interaction of microwave radiation with su- 
perconducting devices has also commanded much inter- 
est. In this paper, we take a slightly different viewpoint 
for these systems by considering them as controllable 
(quantum) optical devices in the microwave range. The 
development of "on-chip superconducting microwave op- 
tical components" (SMOCs) such as phase shifters, mix- 
ers, splitters and detectors, especially in the single and 
few photon regime, will be essential for any future super- 
conducting QIP technology. In such devices nonlinear- 
ity is a fundamental necessity and, despite strong recent 
advances, the resulting quantum dynamics is relatively 
poorly understood and exploited. Our intention in this 
paper is to highlight the potential and the challenges of 
this field through the study of a relatively simple generic 
SMOC device which could, in principle, be developed to 
address the need for a variety of SMOC components. 

Nonlinear effects have recently become a strong focus 
in the area of superconducting quantum circuits because 
of unprecedentedly strong coupling parameters, tunabil- 
ity and engineering flexibility. In the weak nonlinear 
regime the effects are classified into two types: self-Kerr 
effects which give rise to nonlinear phase shifts in a sin- 
gle mode and cross-Kerr nonlinearities which dispersively 
couple two distinct modes. Self-Kerr effects have been 
used to generate squeezed states [I], for parametric am- 
plification [6] and to perform high-fidelity measure- 
ments of quantum bits [MS] , while cross-Kerr effects have 
been discussed in the context of QND measurements and 
intermode coupling [10H12] . Recently, quantum fluctu- 
ations of a pumped nonlinear Kerr-like resonator were 
measured using a qubit as an in-situ probe [T3]. Nonlin- 
earities have also been exploited in other systems, such as 
nanomechanical devices [14j [15] , Bose- Einstein conden- 
sates [ini E] and cold atoms [T51 HH] ■ In this paper we 



explore Kerr effects in a capacitively coupled multimode 
superconducting coplanar resonator which incorporates 
a seven-element SQUID array located at the mid-point 
in the central conductor of the resonator. Figure [T] shows 
optical and electron micrograph images of our sample. 
Different measurements of the same sample can be found 
in [20]. 

Superconducting coplanar resonators (without 
SQUIDs) may be thought as analogous to optical 
Fabry-Perot resonators in terms of mode structure, but 
they exhibit weak Kerr nonlinearities which may cause 
the modes to soften or harden. The strength of Kerr 
nonlinearities in such devices may be greatly increased 
and controlled by inserting nonlinear, nondissipative 
lumped elements based on Josephson junctions [21] . The 
resonator exhibits electromagnetic standing-wave modes 
with a current node at the location of the capacitors. 
Antinodes of current occur at the location of the SQUIDs 
for odd-numbered modes. The instantaneous inductance 
of the SQUID array depends on the total instantaneous 
current through the junctions which may arise from 
both externally applied flux and the presence of photons 
in any of these modes. This dependence induces nonlin- 
earities in the odd-numbered modes as well as couplings 
between these modes. The even-numbered harmonics, 
in contrast, possess nodes of electrical current at the 
location of the SQUIDs and do not exhibit enhanced 
nonlinearities. The SQUIDs also enable flux-tunability 
of the frequencies of the odd harmonic modes, as well as 
enhancement and tunability of the Kerr nonlinearities. 

In order to fully exploit the sensitivity of our SMOC 
we first use the self-Kerr effect to operate one mode as 
a cavity bifurcation amplifier |21j . In this paper we use 
the sensitivity of this mechanism in conjunction with the 
cross-Kerr coupling to demonstrate spectroscopy of other 
cavity modes. It is well known in superconducting qubit 
and SQUID devices that noise sources, thought to be 
embedded in the substrate or the device itself, limits the 
sensitivity and operation of these devices [2"2Tf2~4'] . Elimi- 




FIG. 1: (Color online.) The experimental arrangement for 
pulsed measurements of a nonlinear superconducting res- 
onator intended for microwave optics experiments. RF pulses 
with the envelope shown in insert d are fed to the sample 
(insert a) located in a dilution refrigerator. The power and 
the frequency of the pulses are tuned so that the bifurcation 
regime of the nonlinear resonator is probed. The transmit- 
ted signal is homodyne detected and recorded, insert e shows 
the signal obtained when the sample does or does not exhibit 
bifurcation. Inserts b and c show the input and output ca- 
pacitance and the SQUID array respectively. 



nating these noise sources is one of the grand challenges in 
the field. We observe that our SMOC is also affected by 
a source of low frequency noise and we investigate some 
of its properties. Having characterised the sensitivity of 
our device, we then compare our results to the calcula- 
tion by Dykman |25| , and show that our device operates 
only a factor two worse than this theory predicts. 

Our SMOC consists of a 200 nm sputter-deposited Nb 
film fabricated through optical lithography into a me- 
andering half-wavelength coplanar waveguide resonator 
with end capacitors of design value C — 7pF. The fun- 
damental frequency is ~ 1.77 GHz. An array of N = 7 
Al/AlOx/Al SQUIDs was formed with standard e-beam 
lithography and double-angle shadow evaporation. The 
areal dispersion of the SQUID array is 4%, according to 
SEM imaging. The purpose of using a short array is to be 
able to vary the critical current and hence the magnitude 
of the nonlinearities at the design stage, while maintain- 
ing a fixed overall inductance of the array. The sample is 
magnetically and electrically shielded and attached to the 
mixing chamber of a dilution refrigerator. The measure- 
ment circuit is shown in Fig. [I] Carefully controlled mi- 
crowave pulses are transmitted through the sample using 
lines designed to avoid thermal and amplifier noise from 
reaching the sample. The pulses are homodyne detected 



and measured using an oscilloscope or analogue-to-digital 
converters. We use a pulse repetition rate (5 kHz) that 
is sufficiently slow that the mode relaxes to its ground 
state between pulses. For CW measurements the room 
temperature electronics is replaced by an Anritsu Vec- 
tor Network Analyzer. The measurement bandwidth is 
limited by the circulators to approximately 4 — 8 GHz, 
covering the 2 nd , 3 rd and 4 th modes of the SMOC. The 
third mode is used as our bifurcation detector, the 2 nd 
and 4 th modes are uncoupled and play no further part in 
this discussion. 

We write the Hamiltonian of the system, consisting of 
the 3 rd mode and some other odd n th mode to which it 
couples, as H — H n + H 3 + H : 
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a n (a;) represents the annihilation (creation) operator of 
the n mode of resonant frequency v. t 
Kerr parameter of the n th mode and A 



n.m 
th 



K n is the self- 
is the cross- 



Kerr coefficient between the n th and m th modes. Here 
we assume that all self-Kerr terms except for K% may 
be neglected, as these modes are either operated at low 
power far from bifurcation or are empty. K 3l the self- 
Kerr coefficient of the 3 rd harmonic, is related to the 
critical number of photons iVg required for bifurcation 
via N 3 = 2 r y 3 /y/3K 3 , with 73 being the linewidth of 
the mode. A3.„, the cross-Kerr parameter, defines the 
frequency shift of the 3 rd mode per photon present in 
the n th coupled mode. 

The Kerr coefficients are related to the circuit param- 
eters, and to each other, by 






X 
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P 2 hv ± 
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where Ej is the Josephson energy of a single SQUID of 
inductance Lj and N is the number of SQUIDs in the 
array. /3 — L array / 'L tot is the ratio between the total in- 
ductance of the SQUID array (L array — NLj) and the 
total inductance of the resonator {Ltot = L wg + L array , 
where L wg is the inductance of the waveguide). Both 
self- and cross-Kerr coupling parameters in this device 
are flux-tunable, though not independently. This arises 
from the fact that /?, Ej, v 3 and v n can all be tuned by 
varying the magnetic flux through the SQUIDs. Remark- 
ably, Eq. (4b has an upper limit given by ~ 2-7T J^- « 0.01 
where Zq is the characteristic impedance of the waveg- 
uide and Rk is the resistance quantum. 

To extract some basic parameters of the system, 
we used continuous-wave measurements to acquire the 



forward-scattering parameter, S21, for the 2 nd , 3 rd and 
4 th modes at different values of magnetic flux. As ex- 
pected, the flux modulations of the uncoupled modes 
v-i and 2/4 are negligible, while v 3 shows a modulation 
of ~ 30%. We obtain the following physical parame- 
ters at $/$ = 0: L array = (0.34 ± 0.02) nH, the av- 
erage critical current of a SQUID I c ~ 6.72 /iA, j3 — 
(2.54 ± 0.02)%, v 3 ~ 5.32 GHz and K 3 ~ 940 Hz. Hence 
K 3 /v 3 ~ 2 x 10 -7 , which, given the linewidth of the 
mode in the linear regime 73 ~ 212 kHz, establishes that 
N£, the critical number of photons required for bifurca- 
tion, is ~ 260, which is in reasonable agreement with 
the estimated microwave power at the sample. When a 
cavity mode, in our case the 3 , is probed with a num- 
ber of photons greater than the critical value, the res- 
onator response bifurcates (stochastically and hystereti- 
cally) from the low- to the high-amplitude metastable 
state (see Fig. life) . The probability of the system switch- 
ing (Pb) depends weakly on the quality factor of the 
mode and on the self-Kerr coefficient K 3 while it de- 
pends strongly on the power and frequency of the prob- 
ing pulses, and on the resonant frequency v 3 26J. Fig- 
ure [2] shows 50 bifurcation probability curves as a func- 
tion of the driving frequency. The 10%-90% width of 
the S-shaped averaged probability curve (white line) is 
AS = (4.5 ± 0.5) kHz or AS/v 3 ~ 0.9 ppm. This value 
is in reasonable agreement with the theoretical value cal- 
culated from Dykman's model (AS/v 3 ~ 0.5 ppm) given 
by 
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where 5v = V4 — v 3 is the frequency detuning of the mi- 
crowave driving and T e ff is the effective temperature. 
For comparison, the linewidth of the same mode in the 
linear regime is 73 = 212 kHz or "f 3 /v 3 = 40 ppm. AS* 
represents the sensitivity of the device - i.e. it expresses 
the ability to resolve a frequency shift of AS with a high 
level of confidence within a single shot measurement. 
This bifurcation technique, known as a cavity bifurca- 
tion amplifier [21], is a sensitive threshold measurement 
which may now be used to investigate other aspects of 
the resonator. 

To demonstrate the potential utility of this SMOC, we 
performed spectroscopy of nearby coupled modes. We 
first tuned the bifurcation amplifier to a bifurcation prob- 
ability of 10% and then injected a low-power continuous- 
wave signal and swept its frequency in the vicinity of the 
coupled mode to be detected. As shown in Fig. [3j we 
were able to detect the change in photon occupation of 
the 1 st , 5 th , 7 th and 9 th modes via a variation in the bifur- 
cation probability of the 3 rd mode, thus demonstrating 
the cross-Kerr coupling between the modes. The reso- 
nant frequencies of the coupled modes are obtained via 
a Lorentzian fit providing a value for /3 at zero magnetic 
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FIG. 2: Measurements of the bifurcation probability of the 
3 rd mode of a superconducting nonlinear resonator as a func- 
tion of the driving frequency at 8mK and "I?/$o = 0. Each 
point represents the probability of bifurcation of the mode as 
determined from 1000 pulses. 50 separate consecutively mea- 
sured curves are shown. The white line is the average of the 
50 curves. The expected statistical error of the white line is 
approximately the line thickness. 
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FIG. 3: (Color online.) Spectroscopy of Kerr-coupled mi- 
crowave modes. The bifurcation probability of the 3 rd mode 
is measured while simultaneously applying a continuous swept 
microwave excitation to the 1 st , 5 th , 7 th , and 9 th modes. 



flux of /3 = (2.55 ± 0.1)% that is compatible with the 
previously extracted value. In this way, the cross-Kerr 
effect enables spectroscopy of modes that lie outside our 
measurement bandwidth and thus are undetectable via 
the usual transmission detection. Based on the estimate 
-^- ~ 10, our device is sensitive to ~ 10 photons in the 

1 st mode. The 9 th mode is split for reasons that we do 
not understand; possibly the mode is bifurcating. 

We now consider in more detail the sensitivity of our 
SMOC. We have already noted that Dykman's calcula- 
tion [5S] predicts a width AS due to quantum fluctu- 
ations that is below our experimental value. The bi- 
furcation curves repetitively acquired under nominally 
identical conditions (Fig. [2| show that the bifurcation 
probability Pb is affected by fluctuations that signifi- 
cantly exceed the expected statistical fluctuations. A 



noise source, whose low frequency components appear 
in these measurements, may also be present at higher 
frequencies where it may increase AS and be the source 
of the discrepancy between our experimental values and 
the theory of Dykman. Such a noise source would be 
parametrically coupled to our CJBA, either through the 
frequency of the cavity or through the amplitude and fre- 
quency of the biasing pulses. Estimates of the extrinsic 
sources of noises, such as temperature fluctuations and 
power (frequency) fluctuations of the biasing pulses, sug- 
gest that they cannot account for the observed Pb fluctu- 
ations. Fluctuations in the microwave pulse power at the 
sample holder, measured with the cryostat warm, were 
found to contribute < 10% to the observed AS. Hence, 
in the following, we focus our attention on noise sources 
intrinsic to the device. 

The upper panel in Fig. [4] shows the dependence of 
the width AS* on applied magnetic flux. We first note 
that AS displays a minimum of ASo ~ 4.5 kHz. We 
achieve a reasonable fit by assuming that ASo is an ad- 
ditive constant at all flux values and that the additional 
increase of A S away from zero flux is caused by a quasi- 
static flux noise. The RMS amplitude of the flux noise, 
integrated over the bandwidth of the experiment, can 
be extracted using the measured dependence of 1/3 on 
magnetic flux. This leads to an RMS noise amplitude 
of ~ 5/i$0: which is comparable to commercial SQUID 
amplifiers and probably indicates the quality of our mag- 
netic shielding. The second-order contribution to AS 
from this flux noise at zero flux is ~ 1Hz, which is an in- 
significant contribution to the observed AS'o = 4.5 kHz. 
Hence, AS at zero flux is not due to flux noise. The 
center panel in Fig. [4] shows measurements of AS as a 
function of the refrigerator temperature over the range 
8 — 400 mK at zero magnetic flux. The expected depen- 
dence of bifurcation width with temperature is given by 
Eq|5|, with T eff = |^ coth ^^ (shown as the dotted 
line m the central panel of Fig. El). The effective tem- 
perature T e ff is the physical temperature if T > T co or 
hvs/2kb for T <^i T co , where T co — hv^/Akb ~ 70 mK is 
the crossover temperature from the classical regime to 
the quantum regime [25] . We observe that the values for 
AS are generally within a factor two of Dykman's pre- 
dictions and follow the weak temperature dependence of 
the theory down to the lowest temperature. The lower 
panel in Fig. [4] shows AS as a function of the microwave 
power in units of the critical power for bifurcation. AS is 
found to be almost constant over one order of magnitude 
above the critical power and is again within a factor two 
of the predictions of Dykman's model. 

In summary, we have created a device to assess the 
usefulness of on-chip superconducting microwave optical 
components, or "SMOCs". By inserting a controllable, 
nonlinear element into a superconducting microwave res- 
onator, we are able to create a SMOC that simultane- 
ously exploits the self- and cross-Kerr effects. We have 
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FIG. 4: The width AS of the bifurcation curve as a function 
of reduced magnetic flux (top panel), temperature (central 
panel), and driving microwave power (lower panel). In the 
upper and lower panel the data was taken at refrigerator base 
temperature « 7mK. For the center panel and lower panels 
the data was taken at zero magnetic flux and the dashed lines 
represent the prediction of Dykman's theory. 



shown that high-sensitivity spectroscopy of nearby modes 
may be performed through measurements of the bifurca- 
tion probability of another mode via nonlinear coupling. 
Analysis of the device sensitivity suggests that there is a 
source of noise that manifests as fluctuations in the cavity 
resonant frequency, it has a low frequency component and 
it arises on-chip. The noise may be related to the ubiqui- 
tous noise seen in similar SQUID and qubit systems; we 
have excluded microscopic flux noise and the remaining 
sources include two-level fluctuators in the environment 
and critical current noise in the tunnel barriers. How- 
ever, we note that the intrinsic sensitivity of the device 
is close to the limits given by the theory of Dykman. 
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